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1. Introduction 

Polymers in dilute solution are known to undergo a collapse transition as the 
temperature or solvent quality is changed at what has come to be known as the 0- 
pointfl]. Using universality arguments, it is reasonable to expect the thermodynamic 
behaviour of a lattice model to be the same as the continuum real system as long as the 
dimension of the system, basic symmetries and range of interactions are the same. Such 
lattice models (self-avoiding walks) have been used as models for real, linear polymers 
in solution for over three decades[2]. 

The quality of the solvent may be introduced by the inclusion of short-ranged 
interactions in the model; typically an attractive energy is included for non-consecutive 
nearest-neighbour occupied lattice sites. This model is the standard Interacting Self- 
Avoiding Walk model (ISAW) or O-point modelfl, 3]. The model has been shown 
to accurately predict the critical behaviour of a wide range of real linear polymers in 
solution, not only in the high-temperature phase, but also at the collapse transition, 
which occurs as the temperature is lowered, at the G temperature. The model is 
successful because it captures the strong entropic repulsion between different portions 
of the polymer chain (the self-avoidance), as well as the effect of the difference of 
affinity between monomer-monomer contacts and monomer- solvent contacts (attractive 
interaction). 

Whilst the relevant physical dimension in polymer physics would usually be d = 3, 
the ISAW model has been much studied in two dimensions. This is partly motivated by 
the realisation that d = 3 is the upper critical dimension of the collapse transition, and 
that the model in two dimensions provides an interesting playground. In this paper we 
shall concentrate on the two-dimensional square lattice. 

In the late eighties and early nineties there was much discussion about the 
universality class of the ISAW model, particularly with respect to the adsorption of 
the collapsing walk in the presence of an adsorbing wall [4, 5, 6, 7, 8, 9, 10]. For a while 
there was an apparent contradiction between a slightly modified walk model on the 
hexagonal lattice (the G' model) and the standard G model[4]. This contradiction arose 
in the surface exponents; the exact surface exponents from the G' point model were 
not the same as those calculated numerically for the G modelfll, 12]. The apparent 
contradiction was resolved when it was realised that the exact solution of the G' model 
gives the exponents at the so called special point (where collapse and adsorption occur 
simultaneously) whilst the numerical calculations were at the ordinary point (where 
collapse occurs in the presence of the wall, but without adsorption) [9]. This was verified 
for both models at the two different points using exact enumeration[10]. 

The debate over the G and G' models opened the debate over to what extent the 
nature of the collapse transition depends on the details of the model. Different models 
were examined, and a range of collapse transitions were observed. Blote and Nienhuis 
introduced an 0(n) symmetric model which in the limit n — > gives a bond self- 
avoiding walk model, which is allowed to visit sites more than once but not the lattice 
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bonds[13]. The walk is not allowed to cross. Since the interactions are on the lattice 
vertex, we shall henceforth refer to this model as the vertex-interacting self-avoiding 
walk (VIS AW). This model was shown to have a different collapse transition than the 
point model, with different bulk exponents[14]; the correlation exponents v — 12/23 
for the VIS AW compared to 4/7 for the IS AW and the entropic exponent 7 = 53/46 
compared to 7 = 8/7 for the IS AW. These exponents are conjectured based on a Bethe- 
Ansatz calculation of Izergin-Korepin model[14], and to the best of our knowledge have 
not been numerically tested since their conjecture. 

In recent years there has been a revival in another model with vertex interactions: 
the interacting self-avoiding trails (ISAT) model[15]. This model corresponds to the 
VISAW in which the no-crossing constraint is relaxed. Evidence was presented by one 
of us that the v exponent was also given by v = 12/23, whilst 7 = 22/23(16]. A similar 
situation occurs in the ISAW on the Manhattan lattice, where the walk can only go one 
way down a row or column of lattice bonds, but the allowed direction alternates from 
one row (column) to the next. Here too the correlation length exponent is the same as 
the normal ISAW one, but 7 = 6/7 rather than 8/7[17, 18]. 

Recently the surface exponents of the ISAT model were calculated using transfer 
matrix calculations[19]. We propose here to similarly calculate the surface critical 
behaviour of the VISAW model. In the case of the VISAW model, the no-crosssing 
constraint allows us to extend the transfer matrix calculation using the related density 
matrix renormalisation group (DMRG) method introduced by White[20, 21], applied to 
two-dimensional classical spin models by Nishino[22] and extended to self- avoiding walk 
models by Foster and Pinettes[23, 24]. 

The finite-size calculations rely on results from conformal invariance, which lead one 
naturally to calculate the scaling dimensions x a and x e with fixed boundary conditions. 
Translating these to the more standard exponents requires a knowledge of v. The value 
of v arising from the transfer matrix calculation is at variance with the conjectured 
exact value for the model. We take the opportunity of extending the original transfer 
matrix calculation by Blote and Nienhuis[13]. We find that, up to the lattice widths we 
attain, our best estimate for x e = 0.155 as found by Blote and Nienhuis[13] rather than 
the required x £ = 1/12 = 0.083333 • • •. We conclude that either the finite-size effects are 
particularly severe with this particular model, or a more subtle effect is at play. Either 
way more work is required. 

2. Model and Transfer Matrix Calculation 

The model studied here is defined as follows: consider all random walks on the square 
lattice which do not visit any lattice bond more than once. The walk is not allowed 
to cross at doubly visited sites but may "collide" . A collision is assigned an attractive 
energy —e. The walk is allowed to touch, but not cross, a surface defined as a horizontal 
line on the lattice. Each step along the surface is assigned an attractive energy —e$- 
For the transfer matrix calculation that follows, it is convenient to consider a strip of 
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Figure 1. A vertex interacting self-avoiding walk model showing the vertex collisions, 
weighted with a Boltzmann factor t = cxp(/3e). Surface contacts are weighted 
u s = exp(/3e s ) and a fugacity K is introduced per walk step. The walk is shown 
on a strip of width L = 5. 




width L with an attractive surface both sides of the strip. This is not expected to 
change the behaviour in the thermodynamic limit L — > oo; the bulk critical behaviour 
should not depend on the boundary conditions and when calculating the surface critical 
behaviour, a walk adsorbed to one surface needs an infinite excursion in order to "see" 
the other surface. Additionally, the finite-size scaling results which link the eigenvalues 
of the transfer matrix to the scaling dimensions x s a and x s £ (see (17) and (18)) rely on 
the conformal mapping of the half plane (one adsorbing surface) onto a strip with two 
adsorbing surfaces [25]. A typical configuration is shown in Figure 1. 
The partition function for the model is 

Z= J2 K n u? s t Ni , (1) 

walks 

where K is the fugacity, u s = exp((3es) and r = exp(/3e). N is the length of the walk, 
Ns is the number of steps on the surface, and Nj is the number of doubly-visited sites. 
The average length of the walk is controlled by the fugacity K through 

As K increases from zero, (N) increases, diverging at some value K = K*(u 8 ,t). To 
start we consider what happens in the absence of the adsorbing boundary. For r small 
enough, 

(N)~(K*(u s ,t)-K)-\ (3) 

whilst for large enough r the divergence is discontinuous. Whilst (N) is finite, the 
density of occupied bonds on an infinite lattice is zero, whilst once (N) has diverged, 
the density is in general finite. For small enough r the density becomes non-zero 
continuously at K = K* and for large enough r the density jumps at K — K*. K* may 
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then be understood as the location of a phase transition, critical for r < r co n and first 
order for r > r co ii- The problem of infinite walks on the lattice is equivalent to setting 
K = K* and varying r, then it may be seen that for r < r co n the density is zero and is 
non-zero for r > t co \\. It then follows that t co ii defines the collapse transition point. 

Now let us consider the effect of the adsorbing boundary at constant r. For uj s 
small, the entropic repulsion of the wall is strong enough for the walk to remain in the 
bulk. Once u s is large enough for the energy gain to overcome the entropic repulsion, 
the walk will visit the boundary a macroscopic number of times, and the walk adsorbs to 
the surface. These two behaviours are separated by u s = cj*. For u s < u* the behaviour 
of the walk is not influenced by the wall, and K* is independent of u s . The transition 
K = K* if critical (r < r co n) corresponds to ordinary critical behaviour. However, for 
u s > u*, K* is a function of u s , and the transition is referred to as a surface transition. 
The point K = K*, u s = u* is referred to as the special critical point (again r < r co n). 

As the critical value K* is approached, and in the absence of a surface, the partition 
function (1) and the correlation length £ are expected to diverge, defining the standard 
exponents 7 and v: 

i~\K-KX v (4) 
Z~\K-K*\-^ (5) 

The effect of the surface on the walk is to introduce an entropic repulsion, pushing the 
walk away from the surface. The number of allowed walks is reduced exponentially if 
the walk is constrained to remain near the surface, in particular if one or both ends of 
the walk are obliged to remain in contact with the surface. In this case the divergence 
of Z is modified, and two new exponents are introduced, 71 and 711. Defining Z\ and 
Z\\ as the partition functions for a walk with one end, and both ends, attached to the 
surface respectively, then: 

Z x ~ \K - K*\-^ (6) 
Z11 ~ \K-K*\ ~ 711 (7) 

Whilst the bulk exponents, such as v and 7, are the same at an ordinary critical point 
and at the special critical point, the surface exponents 71 and 711 differ. The exponents 
v, 7, 71 and 711 are related by the Barber relation[26]: 

1/ + 7 = 271 -711. (8) 

The partition function may be calculated exactly on a strip of length L x — > 00 and 
of finite width L by defining a transfer matrix T. If periodic boundary conditions are 
assumed in the x-direction, the partition function for the strip is given by: 

Z L = lim TrfrM. (9) 

L x — >oo V / 

The free energy per lattice site, the density, surface density and correlation length for 
the infinite strip may be calculated from the eigenvalues of the transfer matrix: 

f(K,u 8 ,T) = ^\n{Xo), (10) 
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(11) 
(12) 

(13) 



where Ao and Ai are the largest and second largest (in modulus) eigenvalues. 

Our first task is to find estimates of K*(u s ,t). An estimate for the critical 
point where the length of the walk diverges may be found using phenomenological 
renormalisation group for a pair of lattice widths[28], L and L'. The estimated value of 
K* is given by the solution of the equation: 

k _ | (14) 

Both these methods give finite-size estimates K^(u s ,t) which should converge to 
the same value in the limit L — > oo. Using Equation (4) at the fixed point defined 
by Equation (14), estimates for v and the corresponding surface correlation length 
exponent, u s , may be calculated using 

iu & V dK I dK 



1 
1 



log(L/(L + l)) 

(£/- 



- 1, 



1. 



(15) 



(16) 



u.(L) log(L/(L + l) 
The critical dimensions of the surface magnetic and energy fields may be calculated 
from the first few eigenvalues of the transfer matrix: 

Lin ^ 

(17) 



7T 



Lin 



7T 



with A 2 the eigenvalue with the third largest absolute value. 

The surface scaling dimensions x s a and x s e may be related to the surface correlation 
length exponent v s and the exponent 7711, controlling the decay of the correlation function 
along the surface, through standard relations 

(19) 



1 - 

2xi 



The entropic exponent 7 n is related to 7711 through: 
711 = -m). 



(20) 



(21) 



For a more detailed discussion of the transfer matrix method, and in particular how 
to decompose the matrix, the reader is referred to the article of Blote and Nienhuis [13]. 
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Figure 2. The phase diagram calculated using the Phenomenological Renormalisation 
Group equation. The vertical line is placed at the best estimate of the collapse 
transition, expected to be independent of the surface interaction. (Colour online) 




3. Results 

The finite-size results obtained are, where possible, extrapolated on the one hand using 
the Burlisch and Stoer (BST) extrapolation procedure [29] and on the other hand fitting 
to a three point extrapolation scheme, fitting the following expression for quantity Xl' 

X L = X oa + aL~ b . (22) 

Calculating X^, a and b require three lattice widths. The extrapolated values X^ 
clearly will still depend weakly on L, and the procedure may be repeated, however weak 
parity effects can be seen in their values, often impeding further reasonable extrapolation 
by this method. 

3.1. Phase Diagram 

The phase diagram is shown as a function of u s and r, projected onto the K = K*(t, uj s ) 
plane, in Figure 2. K* is determined using equation (14) using two lattice sizes L, L + l. 
The adsorption line is then fixed by the simultaneous solution of (14) for two sets of 
lattice sizes, L, L + 1 and L + 1,L + 2, so that each line requires three lattice sizes. 
The vertical line is fixed from the best estimate for the bulk collapse transition, here 
Toon = 4.69 [23]. 

In the adsorbed phase, shown in the phase digram in Figure 2, the number of 
contacts with the surface becomes macroscopic, scaling with the length of the walk, and 
the density decays rapidly with the distance from the surface. For the 0-point model 
it has been shown that there is another special line in the phase diagram separating 
the collapsed phase in two: for small enough u s the collapsed walk avoids contacts with 
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Figure 3. The order parameter O = p s /(Lp) plotted as a function of w s for r = 6. 
(Colour online) 




the wall, but for higher values of u s the outer perimeter of the collapsed globule wets 
the surface, defining an attached globule "phase" [30] . To investigate the possibility of 
such a phase, we examine the order parameter for the adsorbed phase (Figure 3) and 
the density of interactions one lattice site out from the wall (there are no interactions 
on the wall, since four occupied bonds must collide) (Figure 4). In the 0-point model 
the presence of such a phase manifests itself by a plateau in the order parameter. Such 
a plateau exists, but starts at or below u s = 1, indicating that the globule is probably 
attached for all attractive wall interaction energies. This is consistent with the plots 
of the normalised density of interactions. Both plots show crossings at a value of u s 
consistent with the adsorption transition. We suggest that the entire phase is "surface- 
attached", and so there is no additional line on the phase diagram shown in Figure 2. 

In Table 1 we locate the collapse transition in strips with fixed walls and u s = 1. 
The collapse transition is determined as follows: solutions to the critical line K^(t) are 
found by using the phenomenological renormalisation group on a pair of lattice widths 
(L and L + 1) and looking for crossings in the estimates for v for consecutive pairs of 
L,L + 1. Since v is different at the collapse point than along r < r co n and r > r co u lines, 
these estimates converge to the correct v for the collapse transition. This gives us the 
following estimates: -fr co ii = 0.3408, r co u = 4.69 and ?7j| rd = 1.77, as well as v CQ \\ = 0.572. 
It is noticeable that its value is much closer to that expected for the 0-point model 
(z/e = 4/7) than the predicted value for this model (z/ ( n =o) = 12/23). We shall see 
that, whilst the estimates for the other quantities of interest are remarkably stable, the 
estimates for u co n seem rather sensitive to how they are calculated. We will return to 
this point later. 
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Figure 4. The order parameter for a possible globule attached phase Ox = pn/(Lp) 
is plotted as a function of ui s for r = 6. The density of interactions one row from the 
surface is used since it is not possible to have collisions on the surface because there 
are only three lattice bonds per surface site. (Colour online) 




Table 1. Location of the collapse value of r with u> s = 1, using three lattice widths. 
Estimates for -ftf co ii and v co i\ are also shown, as well as estimates for rj°* d calculated for 
uj s = 1 The three point extrapolations are shown in the second half of the table. 



L 




Tcoll 


^coll 




3 


0.359410 


4.071423 


0.614465 


1.024334 


4 


0.351725 


4.410526 


0.596955 


1.147824 


5 


0.347865 


4.540658 


0.588407 


1.233790 


6 


0.345694 


4.598914 


0.583712 


1.297254 


7 


0.344369 


4.628215 


0.580898 


1.346184 


8 


0.343508 


4.644460 


0.579079 


1.385168 


9 


0.342920 


4.654221 


0.577827 


1.417020 


10 


0.342502 


4.660572 


0.576905 


1.443585 


BST 


oo 0.3408 


4.673 


0.574 


1.77 




Three point extrapolated 


results 




3 


0.339412 


4.696857 


0.571064 


2.103158 


4 


0.340217 


4.681105 


0.572693 


1.963018 


5 


0.340540 


4.676727 


0.573260 


1.901623 


6 


0.340676 


4.676316 


0.573231 


1.868539 


7 


0.340749 


4.676871 


0.573019 


1.845418 


8 


0.340767 


4.678731 


0.572356 


1.832298 


BST 


oo 0.3408 


4.69 


0.572 


1.77 
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Table 2. Special point location with r unconstrained, using four lattice widths. 
Estimates for 77jj p are also shown. The three point extrapolations arc shown in the 
second half of the table. 



L 


K co n 


Tcoll 




sp 


3 


0.335871 


4.989134 


2.452162 


-0.120915 


4 


0.337720 


4.868882 


2.418298 


-0.110883 


5 


0.338679 


4.809216 


2.399537 


-0.104452 


G 


0.339256 


4.774456 


2.387539 


-0.099806 


7 


0.339624 


4.752731 


2.379405 


-0.096306 


Q 
O 


0.339874 


4.738301 


z.o ( ooyo 


-u.uyoooo 


9 


0.340050 


4.728276 


2.369291 


-0.091352 


BST 


oo 0.3408 


4.6901 


2.3513 


-0.07843 




Three point extrapolated 


results 




3 


0.340975 


4.682860 


2.344211 


-0.06881 


4 


0.341078 


4.676046 


2.338622 


-0.059977 


5 


0.340899 


4.684292 


2.343044 


-0.064925 


6 


0.340848 


4.686507 


2.343975 


-0.065385 


7 


0.340811 


4.688205 


2.344840 


-0.066177 



In Table 2 we seek the special point along the collapse transition, in other words 
the point at which the extended, collapsed and adsorbed phases co-exist. A different 
set of critical exponents is expected. In order to find the special point, we need an extra 
lattice width; three are required to find oj*(t) and a fourth is required to fix the collapse 
transition. We find K co n = 0.3408 and r co n = 4.69 as in the case when u s = 1. The 
special point is found to be at a value of = 2.35 ± 0.01. The estimate of ?7jj p is not 
very precise, but here seems to be around —0.06 — > —0.07. 

3.2. Results for the semi-flexible VIS AW 

The model may be extended by introducing different weighting for the corners and the 
straight sections. We follow the definitions in Reference [13] and add a weight p for each 
site where the walk does not take a corner (i.e. for the straight sections). As p is varied 
we expect the collapse transition point to extend into a line. It turns out that there 
is an exactly known point along this line. The location of this point is given exactly 
as[13]: 

z = K 2 t = {2- [l-2sin(0/2)] [l + 2sin(0/2)] 2 }~ 1 ' 

K = -4zsin(#/2)cos(7r/4-fl/4) > ^ 

pK = z [1 + 2 sin(0/2)] 
6 = -tt/4 

This gives exactly the location of the multicritical collapse point when p = p* = 
0.275899 ■ • • as K co]1 = 0.446933 • • • and r coU = 2.630986 • • •. Using this exactly known 
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Table 3. Table shows: 

A the location of the value of r co n with w s = 1 and p — p* , using three lattice widths 
and crossings of v estimates. Estimates for K co \\ and ?yj| rd are also shown. 
B Estimates for oj s for the ordinary point along with the estimates for r)°^ d at this 
point with K = K co n, p = p* and r = r co u. Only two lattice sizes are required to fix 
uj s , hence we can calculate up to a lattice width of L = 11. 



A 



B 



L 


K co \\ 


Tcoll 


?y° rd 


w ord 


77° rd 


3 


0.464018 


2.309912 


1.401892 


0.760808 


1.813498 


4 


0.457207 


2.451700 


1.471983 


0.785333 


1.787052 


5 


0.453616 


2.520291 


1.520159 


0.797646 


1.776227 


6 


0.451604 


2.556015 


1.554337 


0.805442 


1.770439 


7 


0.450391 


2.576330 


1.579631 


0.811096 


1.766807 


8 


0.449615 


2.588789 


1.599034 


0.815550 


1.764286 


9 


0.449089 


2.596962 


1.614375 


0.819245 


1.762416 


10 


0.448720 


2.602583 


1.626774 


0.822417 


1.760965 


11 








0.825204 


1.759802 


BST oo 


0.4473 


2.597 


1.708 


0.8955 


1.7499 


exact 


0.446933 


• • • 2.630986 • • 








Three point extrapolated results 


3 


0.444582 


2.656251 


1.953742 


0.824571 


1.761506 


4 


0.446572 


2.628298 


1.807661 


0.789157 


1.757867 


5 


0.447052 


2.622897 


1.774328 


0.799994 


1.755193 


6 


0.447197 


2.621984 


1.762266 


0.863901 


1.753208 


7 


0.447252 


2.622163 


1.754657 


0.881961 


1.751795 


8 


0.447224 


2.623411 


1.753850 


0.816627 


1.750839 


9 








0.820168 


1.750223 



point we hope to be able to extend the number of different data points and improve the 
precision of the determination for different surface exponents. 

In Table 3 we calculate estimates for ?7j| rd in two ways. Firstly we fix u s = 1 and 
p = p* and determine K* by solving Equation (14) and determining the multicritical 
point by looking for crossings in the estimates for v. Fixing the multicritical point 
this way requires three lattice sizes, L, L + 1 and L + 2. Estimates for K co \\ and t co ii 
calculated in this way are shown in the columns marked A and are seen to converge 
nicely to the expected values. The second method used consisted in fixing K, t and 
p to their exactly known multicritical values, and fixing u s to the ordinary fixed point 
looking for solutions to Equation (14). This only requires two lattice widths, giving an 
extra lattice size. The values of r]^ rd are shown as calculated from the two methods, and 
converge to values consistent with the p — 1 case. 

In Table 4 we present results calculated fixing r = r co u and looking for simultaneous 
solutions of the phenomeno logical renormalisation group equation (14). These solutions 
exist at two values of u s , the ordinary and the special fixed points. The values of ifcoib 
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Table 4. Location and estimates of 7711 for both the ordinary and special point fixing 
t = T co n(p*), and using three lattice widths. The three point extrapolations are shown 
in the second half of the table. 



Ordinary Point 



Special Point 



J 

Li 


-ft coll 


, ,ord 


~iOrd 

"ll 


A coll 


u s 


1™ 


3 


0.444849 


0.727730 


1.876811 


0.444289 


3.840487 


-0.254004 


4 


0.446261 


0.765809 


1.816907 


0.445726 


3.660264 


-0.191557 


5 


0.446626 


0.782852 


1.795039 


0.446279 


3.575039 


-0.157742 


6 


0.446763 


0.792806 


1.784181 


0.446537 


3.527515 


-0.136796 


7 


0.446826 


0.799587 


1.777736 


0.446675 


3.498075 


-0.122630 


8 


0.446861 


0.804688 


1.773440 


0.446754 


3.478454 


-0.112441 


9 


0.446881 


0.808784 


1.770341 


0.446804 


3.464654 


-0.104773 


10 


0.446894 


0.812223 


1.767979 


0.446836 


3.454537 


-0.098797 


BST 00 


0.446933 


0.8529 


1.75 


0.44693 


3.4029 


-0.05241 


exact 


0.446933 • ■ 












Three point extrapolated results 


3 


0.445398 


0.740647 


1.855322 


0.444797 


3.415470 


-0.065260 


4 


0.446899 


0.821222 


1.809568 


0.445915 


3.410716 


-0.062347 


5 


0.446915 


0.830707 


1.791541 


0.446934 


3.506065 


-0.060003 


6 


0.446923 


0.794846 


1.782160 


0.446933 


3.407005 


-0.058240 


7 


0.446928 


0.801100 


1.755376 


0.446932 


3.406126 


-0.056866 


8 


0.446931 


0.867939 


1.753498 


0.446932 


3.405498 


-0.055762 



u s and 7711 are given for the two fixed points. Again agreement is found with previous 
values calculated. 

At the ordinary point the exponent v s is expected trivially to take the value —1, and 
this was verified in the various calculations at the ordinary point, with good convergence. 
At the special point the correlation length along the surface is not expected to be trivial. 
In order to obtain the best estimate for this exponent we determined the location of the 
special point by fixing K, r and p to their multicritical values and then determining 
the special point by looking again for solutions to the phenomenological renormalisation 
equation (14), shown in Table 5 and using Equation (15) and Equation (16). One may 
also obtain an independent estimate for v s calculating the scaling dimension x\ using 
(18), from which v s = (1 — x|) _1 . The special point was determined using the odd sector 
of the transfer matrix, whereas the calculation of x s e requires the even sector, so whilst 
the determination method only gives one estimate of ?7jj p it gives two estimates (one for 
each lattice size) for the critical dimension x s £ . These different estimates are shown in 
Table 5. The values of u s s p converge to 1.487 whilst x £ = 0.332, which gives uf = 1.497. 
It seems likely that uf = 1.49 ± 0.01. Again the estimates of v do not converge to 
v = 12/23, but neither do they converge to the values found above for u s = 1 and p = 1 
(see Table 1). The crossover exponent is calculated from the estimates found for each 
size 4> s = u/u s , therefore the extrapolated value is only as good as the estimated values 
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Table 5. Location of the special point for K = K co u, p = p* and r = r co n using 
intersections of rf^ . The three point extrapolations are shown in the second half of the 
table. The values of x s e are calculated in the even sector at the value of ui s s p . 



L 




"IT 


;/ 




(j) s = v j v s 




xl{L + 1) 


3 


3.575571 


-0.170489 


0.527116 


1.829877 


0.288061 


0.527910 


0.422737 


4 


3.506382 


-0.135551 


0.534983 


1.725710 


0.310008 


0.449880 


0.401500 


5 


3.472677 


-0.116031 


0.540007 


1.668509 


0.323646 


0.416568 


0.388462 


6 


3.453634 


-0.103697 


0.543502 


1.632771 


0.332874 


0.397954 


0.379518 


7 


3.441756 


-0.095233 


0.546069 


1.608442 


0.339497 


0.386019 


0.372974 


8 


3.433803 


-0.089075 


0.548029 


1.590960 


0.344477 


0.377697 


0.367970 


9 


3.428190 


-0.084396 


0.549571 


1.577675 


0.348349 


0.371554 


0.364018 


10 


3.424062 


-0.080719 


0.550811 


1.567260 


0.351393 


0.366832 


0.360818 


11 


3.420926 


-0.077751 


0.551828 


1.558893 


0.353944 


0.363087 


0.358175 


BST 


oo 3.402 


-0.0499 


0.5592 


1.487 


0.3777 


0.332 


0.332 




Three point extrapolated 


results 










3 


3.404662 


-0.056758 


0.567274 


1.504747 




0.367160 


0.330606 


4 


3.404891 


-0.056025 


0.566011 


1.501149 




0.351913 


0.326878 


5 


3.404606 


-0.054962 


0.564760 


1.497315 




0.344469 


0.326617 


6 


3.404268 


-0.053981 


0.563832 


1.498468 




0.340278 


0.327150 


7 


3.403967 


-0.053136 


0.563148 


1.488934 




0.337719 


0.327793 


8 


3.403715 


-0.052418 


0.562652 


1.488636 




0.336068 


0.328375 


9 


3.403508 


-0.051807 


0.562224 


1.488766 




0.334959 


0.328865 



of v and v s . If we believe v = 12/23 and v s = 1.5 then (f> s = 8/23 = 0.34782 ■ ■ -. 
4. Extending the results with DMRG 

One of the limitations of the transfer matrix method is the limited number of lattice 
widths that may be investigated. One way of getting round this problem is to generate 
approximate transfer matrices for larger widths. There exists a method of choice for 
doing this; the Density Matrix Renormalisation Group Method (DMRG) introduced by 
White[20, 21], extended to classical 2d models by Nishino[22] and self-avoiding walk 
models by Foster and Pinettes[23, 24]. 

The DMRG method constructs approximate transfer matrices for size L from a 
transfer matrix approximation for a lattice of size L — 2 by adding two rows in the middle 
of the system. This process is clearly local, whereas the VISAW walk configurations are 
clearly non-local. This problem is solved by looking at the model as the limit n — > of 
the 0(n) model. 

The partition function of the 0(n) model is given by: 

Z 0{n) = Y.n l K N p N ^T N \ (24) 
Q 

where Q denotes the sum over all graphs containing loops which visit lattice bonds at 
most once and which do not cross at lattice sites and / is the number of such loops. N st 
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Figure 5. Local complex vertices for the DMRG method 
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is the number of straight sections. In the limit n — > the model maps onto the expected 
model with the odd sector of the corresponding transfer matrix giving the walk graphs, 
as above. Viewing the model in this way enables us to map the loop graphs into oriented 
loop graphs. Each loop graph corresponds to 2 7Vlo °p s oriented loop graphs. We associate 
different weights n + and n_ for the different orientations, with n = n + + n_, and this 
enables us to rewrite the partition function as follows: 

Zo(n) = E (*+ + ™J K n p n ^t n i (25) 
g 

= £ n l +n l I K N p N *u? s T Nl . (26) 
Q* 

Whilst the weights n + and n_ are still not local, they may be made local by realising 
that four more corners in one direction than the other are required to close a loop 
on the square lattice. If we associate a with each clockwise corner and a~ l for each 
anti-clockwise corner we find n + = a 4 and n_ = a -4 , setting a = exp(i9/4) gives: 

n = a 4 + a~ 4 = 2cos#. (27) 

The model studied here then corresponds to 9 = n/2. The resulting local (complex) 
weights are shown in Figure 5. 

Now that the vertices are local, the DMRG method may be applied. The vertices 
represented in Figure 5 may be most easily encoded by defining a three state spin on 
the lattice bonds, for the horizontal bonds the three states would be arrow left, empty 
and arrow right. For details of the DMRG method the reader is referred to the articles 
by White[20, 21] and Nishino[22], but in essence the method consists in representing the 
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Figure 6. Figure shows the schematic transfer matrix obtained from DMRG iteration. 
Circles show spins defined for the original model (3 state spins for the lattice bond: 
empty, and two arrow states). Squares show the m-state pseudo spins. The open 
circles are summed. The projection of the upper half is also shown schematically. 



m 




transfer matrix for the VISAW model as the transfer matrix for an equivalent system 
where the top and bottom of the strip is represented by an m-state pseudo-spin with 
only the inner two of rows kept explicitly in terms of the original 3 state spins. For 
small lattice widths this identification may be done exactly, and the interaction matrix 
may be chosen exactly, but for a fixed value of m there will come a stage where this 
procedure is no longer exact. At this stage the phase space in the m-spin representation 
is smaller than for the real system and an approximation must be made. Starting from 
the largest lattice width that may be treated exactly by the pseudo-spin system, two 
vertices are inserted in the middle (see Figure 6). The 3 x m states at the top of the 
system must be projected onto m states to recover a new pseudo-spin system. This must 
be done so as to lose the smallest amount of information, and this is where the DMRG 
method comes in. It turns out that the best change of basis is derived by constructing 
the density matrix for the top half of the lattice strip from the ground-state eigenvector 
by tracing over the lower half system. The density matrix is then diagonalised and the 
m basis vectors corresponding to the m largest eigenvalues of the density matrix are 
kept. 

There are two modifications on the basic method which improve the quality of 
the results obtained. The number of left arrow minus the number of right arrows is 
conserved from one column to the next, so the transfer matrix may be split into sectors 
which are much smaller than the original. Since the DMRG method may be viewed as 
a variational method, the quality of the results may be improved by using the scanning 
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Figure 7. Location of the ordinary collapse point (A) and corresponding value of r]^ rd 
(B) with p = p* , K = K co \\ and r = r co n. (Colour online) 
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(or finite size method) DMRG method where once the desired lattice width is obtained 
one grows one half of the system and shrinks the other, whilst projecting as before, so 
that the exactly treated spins move across the system. As few as three or four sweeps 
is known to vastly improve the precision of the method[20, 21]. 

Clearly the precision of the method is controlled by m; the larger m the greater 
the information kept. In what follows we varied m up to values of m = 200 and verified 
that good convergence was obtained. This conditioned the lattice size we looked at in 
DMRG. Whilst physical quantities such as the density converge rapidly with m, scaling 
dimensions (which interest us here) converge more slowly. As a result the largest lattice 
width presented here is L = 20, which nevertheless corresponds to a good improvement 
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Figure 8. Location of the special collapse point (A) and corresponding value of rf^ 
(B) with p = p* , K = K co \\ and r = r co ii. (Colour online) 
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over the pure transfer matrix method. 

For the DMRG calculation we fixed p = p*; K = K co n; r = t co \\ and used the 
solutions of Equation (14) to find the ordinary and special fixed points as well as the 
corresponding r]\\. The x s e were calculated from the even sector at these fixed points. 

In Figure 7 and Figure 8 we show the DMRG results along with the transfer matrix 
results for ujg and r/y for the ordinary and special points. We deduce for the ordinary 
point u;£ rd = 0.86±0.01 and r/f[ rd = 1.75±0.01 and for the special point oof = 3.41 ±0.01 
and r/j p = -0.05 ±0.01. 

In Figure 9 we show the estimates for the scaling dimension x s £ calculated at the 
special point. We determine x s e = 0.333 ± 0.001. This leads to = 1.5. 
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Figure 9. Calculation of x s e using DMRG at the special collapse transition with 
p = p* , K = K co \\ and r = r co ii. (Colour online) 
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Figure 10. Bulk scaling dimensions x a and x e calculated at p = p* , K = K co \\ 
and r = T C oii for periodic boundary conditions using DMRG with m up to m = 190. 
(Colour online) 




5. Discussion 

To conclude, we summarise the exponent values found: 



Method 


ord 








TM 


1.75 ±0.05 


-0.05 -0.08 


1.48 ±0.04 


0.332 ±0.005 


DMRG 


1.75 ±0.01 


-0.05 ±0.01 




0.333 ±0.001 
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As discussed above, the calculation of the exponents 71 and 711 as well as <f> s require a 
knowledge of the bulk exponents v and 7. Whilst there are conjectured exact values 
for these exponents, and in particular the exponent v = 12/23 is found to a good level 
of precision for the Trails model which tends to lend support to this value, the transfer 
matrix calculations for the VISAW walk model do not seem to reproduce the required 
values, Blote and Nienhuis[13] find x e = 0.155 rather than the required x e = 1/12 for 
example. We extend the transfer matrix results for the periodic boundary conditions 
using DMRG, and find a result for x £ consistent with Blote and Nienhuis[13](see 
Figure 10). 

Further work is required to calculate the exponents by different methods, for 
example Monte Carlo, in order to understand the apparent differences in results which 
arise in the exponent v. Either the differences are a result of particularly strong finite- 
size scaling effects, which is surprising since the surface exponents themselves seem to be 
remarkably stable in comparison, or perhaps an indication that the critical behaviour of 
this model is more subtle than initially thought. Either way the model warrants further 
study. 
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